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ABSTRACT 



Context. Direct detection of exoplanets requires high dynamic range imaging. Coronagraphs could be the solution, but their perfor- 
mance in space is limited by wavefront errors (manufacturing errors on optics, temperature variations, etc.), which create quasi-static 
stellar speckles in the final image. 

Aims. Several solutions have been suggested for tackling this speckle noise. Differential imaging techniques substract a reference 
image to the coronagraphic residue in a post-processing imaging. Other techniques attempt to actively correct wavefront errors using 
a deformable mirror. In that case, wavefront aberrations have to be measured in the science image to extremely high accuracy. 
Methods. We propose the self-coherent camera sequentially used as a focal-plane wavefront sensor for active correction and differ- 
ential imaging. For both uses, stellar speckles are spatially encoded in the science image so that differential aberrations are strongly 
minimized. The encoding is based on the principle of light incoherence between the hosting star and its environment. 
Results. In this paper, we first discuss one intrinsic limitation of deformable mirrors. Then, several parameters of the self-coherent 
camera are studied in detail. We also propose an easy and robust design to associate the self-coherent camera with a coronagraph that 
uses a Lyot stop. Finally, we discuss the case of the association with a four-quadrant phase mask and numerically demonstrate that 
such a device enables detection of Earth-like planets under realistic conditions. 

Conclusions. The parametric study of the technique lets us believe it can be implemented quite easily in future instruments dedicated 
to direct imaging of exoplanets. 



Key words, instrumentation: adaptive optics — instrumentation: high angular resolution 
techniques: high angular resolution — techniques: image processing 



instrumentation: interferometers 



1. Introduction 

Exoplanets are typically 10 7 to 10 10 fainter than their host and 
are often located within a fraction of an arcsecond from their star. 
Numerous coronagraphs have been proposed to reduce the over- 
whelming light of a star to obtain a direct imaging of extraso- 
lar pl anets dRouan et all 120001: iMawet et all 120051: iGuvon et all 
2005). Several of them prov ide observations (ISchneider et all 
1998; Bo ccaletti et all |2004|). But performance is limited by 
wavefront errors in the upstream beam for all these corona- 
graphs and the final focal plane image shows stellar speckles. 
The effect of most of these aberrations ca n be corrected by adap- 
tive optics (AO) or eXtreme AO (XAO. IVerinaud et all [2008) 
but the uncorrected part generates quasi-static speckles, which 
limit the image contrast dCavarroc et all l2fJ06l: [Macintosh et all 
120051) . To reduce this speckle noise, differential imaging tech- 
niques attempt to subtract a reference image of the stellar speck- 
les from the science image (star plus companion). 

Several ways ar e used to m e asure t hat reference : spec- 
tral c haracteristics dRacine et al. , 19991: Marois et all 2000. 
20041) . polarization states ([Baba&Murakami, 2003; S tam et all 
20041). differential rotation in image ( ISchneider et al [1998; 



Maroi s et a [2006), or incoherence between stellar and com- 
panion lights dGuvoni |2004|) . The self-coherent camera with 



which we work uses the last property {Baudoz et al., 2006; 
Galicher & Baudoz, 2007). But before using one of these a pos- 
teriori techniques, we may actively correct quasi-static wave- 
front errors so that a first speckle reduction is achieved and dif- 
ferential imaging techniques have less to do. Because of the 
low level of aberrations that must be achieved (a few nanome- 
ters), the wavefront sensor of that loop has to measure for phase 
and amplitude errors in the final science image to avoid differ- 
ential errors introduce d by classical wavefront sensor (Shack- 
Hartm ann for example, IShack & Platl 1197 ll) . ICodona & Angell 
d2004l) suggest using the incoherence between companion and 
stellar lights and use a modified Mach-Zender interferometer to 
encoded the stellar speckles. 

The instrument we propose, the self coherent camera (SCC), 
is based on the same property and uses Fizeau interferences. We 
insist on how SCC can be used both as a wavefront sensor for 
active correction (called step A in thi s paper) and as a d iffer- 
ential imaging technique (step B). In Galic her et al.l d2008l) . we 
numerically demonstrated that, applying step B after step A, a 
self-coherent camera associated with a 32x32 deformable mir- 
ror and a perfect coronagraph detects earths (contrast of 2 10 -10 ) 
from space in a few hours under realistic assumptions (zodiacal 
light, photon noise, read-out noise, phase errors of 20 nm rms, 
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20% bandwidth and 1 % effective bandwidth, 8 m telescope with 
a 25% throughput and Aq = 800 nm). Here, we report results 
from a parametric study of the SCC, and propose an easy and ro- 
bust design to associate it with coronagraphs that use a Lyot stop. 
We explain the performanc e in the case of a four-quadrant phase 
mask coronagraph (FOPM. lRouan et aTll2000h . Section|2]recalls 
the principle of the technique and presents the estimators of the 
pupil complex amplitude (phase and amplitude errors, step A) 
and companion images (step B). Section[3]provides the assump- 
tions and criterions used for the parametric studies. Section [4] 
is a general study (no SCC) of one intrinsic limitation for de- 
formable mirrors and presents the best contrast they can pro- 
vide. The signal-to-noise ratios on both SCC estimators (wave- 
front and companion) are developed in Sect. [5] Section [6] es- 
timates the required stability of the reference beam. We report 
the effects of amplitude aberrations on the SCC performance in 
Sect. [7] The last two sections are the most important in the pa- 
per. The first one is the study of the chromatism on the SCC per- 
formance when a perfect coronagraph is used (Sect. [8]). In the 
second one (Sect. [5J, we propose a device to associate the SCC 
and any coronagraph using a Lyot stop. We discuss the case of 
a FQPM coronagraph: earths are detected using such a corona- 
graph with the SCC just by adding a small hole to the Lyot stop. 

2. Self-coherent camera principle 

The goal of the self-coherent camera is the measurement (step A) 
of phase and amplitude aberrations in the pupil upstream the 
coronagraph and the speckle field estimation (step B) without in- 
troducing any non-common path errors. To do so, we use spatial 
interferences in the science image to encode the stellar speckles 
that are directly linked to the wavefront aberrations. This section 
briefly recalls how to use the self-coherent camera to measure 
wavefront errors and also to reduce the speckle noise in the im- 
age. 

More details c an be found in iBaudoz et alj ((2006); 
Galic her et al] (2008). As in these previous papers, we consider 
hereafter only space observations. Figure Q] presents the instru- 
ment schematics. A deformable mirror, located in a plane con- 



jugated to the entrance pupil, reflects the beam coming in from 
the telescope. We then split the beam. In the image channel, the 
beam goes through a coronagraph. In the reference channel, we 
suppress all the companion light and extract a beam containing 
only light from the host star (cf. Sect. [9j. Finally, we recombine 
the two beams in a Fizeau scheme to obtain spatial fringes in 
the science image on the detector. Phase and amplitude aberra- 
tions give a coronagraphic residue in the last focal plane. And 
the reference channel induces spatial interferences on these stel- 
lar speckles (cf. Fig 13, whereas it does not have any impact on 
a possible companion image since companion light is not coher- 
ent with star light. The stellar speckles are thus encoded (mod- 
ulated), whereas the companion image is not. We propose the 




Fig. 2. Science image of the self-coherent camera. Stellar speck- 
les of the coronagraphic residue are fringed, hence spatially en- 
coded. No companion is present. 
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Fig. 1. Self-coherent camera principle schematics. 



following protocol for detecting faint companions: 

- Step A, wavefront estimation and correction: estimate 
phase and amplitude errors from the focal plane image and 
correct for them using a deformable mirror (correction loop 
without any differential errors, Sect. |2.2| > 

- Step B, companion detection: record the science image 
when the best correction is achieved and post-process that 
image to overcome the DM limitation (Sect. [23T l. 



2.1. Notations 

Before giving the expression of the wavefront estimator, we 
summarize the notations of this paper in Table Q] Their signif- 
icance is given in the text when necessary. 



2.2. Step A: correction loop 

2.2.1 . Estimation of the aberrated wavefront 



In lGalicher et all d2008l) . we proposed an estimator of phase er- 
rors in the pupil upstream of the coronagraph. Here, we give 
the SCC estimator of both phase and amplitude aberrations, in 
other words, the pupil complex amplitude. 
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Table 1. Table of the notations. Index i represents R, S, or C re- 
ferring to the reference channel, the stellar residue of the image 
channel and the companion light of the image channel. 



That amplitude is extracted from the science image / (Fig.|2j. 
Using notations from TableQ] we can write the interferential im- 
age I as 



\A s (a)\ 2 + |A R (a)| 2 + |A c (a)| 2 



+2Re(A s (a)A R (a) exp(^f^)) 



dX, 



(1) 



where the wavelength X belongs to = [Xq - AX/2, Xq + AX/2]. 
The first three terms are the coronagraphic residue (stellar speck- 
les to be reduced), the reference image, and the companion im- 
age (to be detected). They constitute the unmodulated part of the 
interferential image. The last term traduces the spatial speckle 
modulation: interferences with a X/tjo interfringe. The required 
pupil separation (Fig. [TJ not to overlap the peaks in the spatial 
frequency plane (inverse Fourier transform of the image 7) is at 
least: 



ft = 1.1 (1.5 + 0.5/y) Z>. 



(2) 



We use a 10% margin (factor 1.1). 

We follow iGalicher et al.l d2008l) to estimate wavefront er- 
rors, and we extract the modulated part I-(a) of 1(a) because it 
contains a linear combination of As (what we look for) and Ar. 
To do it, we select one of the two lateral peaks of the inverse 
Fourier transform of I and apply a Fourier transform to it: 



I-(a) 



j R l 2 A s (a)Al(a)zxp( 2 J^\n. 



(3) 



We then assume a small effective bandwidth (cf. Sect. [8j and 
estimate the pupil complex amplitude downstream the corona- 
graph: 



<As(«) - T- 



Ua) 



A* R (a)F(a) 



(ii), 



(4) 



where T~ l is the inverse Fourier transform, u the cooordinate 
in the correlation plane of 7, and F is applied to partially com- 
pensate for chromatic effects (see Sect. 18.21 ). Assuming small 
phase (</>) and amplitude (a) aberrations and a perfect corona- 
graph, we find the complex amplitude t//s downst ream of the 
coronagraph (see Eqs. 6 to 8 in lGalicher et alll200 8. for details): 



i//s(u) — tfiQ [a(u) + i (f>(u)] P(u), 



(5) 



where psio is the star flux upstream of the coronagraph. From 
the two last equations, we deduce 



<p(u) — ia(u) =i —if 



I-(a) 



L// A R (a) F(a) 



(a). 



(6) 



We can thus estimate wavefront errors upstream of the corona- 
graph from /_, extracted from the science image. The other terms 
can be estimated as follows. The chromatic factor is computed, 
and its expression is given in Sect. 18.21 The constant i//o is uni- 
form (TableQ. We can precisely estimate it from flux calibration 
of the space telescope. Finally, the reference amplitude in the fo- 
cal plane A R has to be recorded before the beginning of the loop 
and has to be stable. We point out that such a stabilization can 
be achieved (Sect. [6}. 

Once wavefront aberrations are estimated, we correct for them 
with a deformable mirror as explained in the next section. 

2.2.2. Correction of the aberrated wavefront 

In the whole paper, we assume that we use only one deformable 
mirror with N dC oi-N- dC t actuators to correct for the wavefront er- 
rors. If only phase aberrations exist, it is possible to clean all 
the ([-/Vact/2, Afact/2] (Xq/D)) 2 centered area in the science im- 
age. If both phase and amplitude aberrations are present, we can 
only clea n half of that r e gion, and we adopt the method pro- 
posed by Borde & TraubJ d2006l) to make the phase error screen 
hermitian. As the estimation is not perfect (small aberration lin- 
earization, noises, reference division) and the speckles are not 
static but quasi-static, we need to iterate a few times. 

To achieve high-contrast imaging, we need to measure wave- 
front aberrations with high accuracy and to drive the deformable 
mirror faster than the aberration evolution time. As the required 
time to work out the SCC estimator is relatively short (three fast 
Fourier transform), the correction speed is determined by the 
integration time needed to reach a reasonable signal-to-noise 
ratio for the wavefront measure. There is then a compromise 
between the achievable image contrast and the pointed stellar 
flux. That compromise is roughly the same regardless of the 
speckle calibration technique. Recent industrial studies for a 
space-based 1.5m telescope est imate that st abilities can be ex- 
pected as low as 1 pm per hour dGuvonj [20091) . 



2.3. Step B: companion estimation 

Once the best correction is achieved (last iteration of step A), we 
apply the p ost-proccessing algor i thm d escribed in Baudoz et al.l 
(2006) and Galicher & Baudoz (2007) to suppress the largest 
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part of the speckle residue. We estimate the companion im- 
age I c est using 

i r \ i < \ i t \ l 7 -( a )| 2 /n \ 

/CestW = /cem(a) - Ma) - - , (7) 

Ma) 

where /_ is defined by Eq.[3]and / cent is the unmodulated part of 
the science image / (Eq.Q]), which can be written as 



War) = j j2 



h(a) + Ma) + I c (a) 



dX. 



(8) 



To extract it, we select the central peak of the inverse Fourier 
transform of the science image (Eq.Q3 and apply a Fourier trans- 
form to that selection. Finally, /r = |Ar| 2 is the reference image. 
It has to be estimated as explained in the introduction of Sect. [6] 

3. Assumptions and criteria 

This section introduces the assumptions of our numerical studies 
and the criteria for optimizing all the parameters of the instru- 
ment. 



3.1. Assumptions 

In the whole paper, we assume spatial observations (no dynamic 
aberrations) and achromatic coronagraphs without defects (per- 
fect coronagraph or FQPM). The power spectral density (PSD) 
of static phase errors in the instrument upstream of the coron- 
agraph varies as / 3 , where / is the spatial f requency, which 
corre sponds to typical VLT mirror aberrations (Borde & Traub, 
2006). Amplitude aberration PSD is flat or evolves as f~ 3 . The 
reference complex amplitude A R in the detector plane is the 
Fourier transform of 

- the pupil upstream of the coronagraph (same wavefront er- 
rors) densified to obtain a pupil diameter Dr instead of D, 
if a perfect coronagraph is used. The flux is set to verify the 
condition of Sect. l5.ll in the corrected area. 

- a pupil of diameter Dr extracted from the Lyot stop 
plane (phase and amplitude errors depend on the corona- 
graph) if a FQPM coronagraph is used (see Sect. [9]). 

We consider the deformable mirror as a continous face sheet 
supported by N act xN act actuators arrang ed in a square pattern of 
constant spacing dBorde & Tr aub. 2006) and located in a plane 
conjugate to the entrance pupil. The nth-actuator influence func- 
tion is exp (In (0. 15) (N act (f — $ n )/D ) 2 ), where £ n is the center of 
the nth-actuator (Huang et al.l 12008b . We use th e method of en- 
ergy m inimization in the pupil plane (Eq. A3 of Borde & Traub, 
2006) to project the estimated phase on the deformable mirror 
and shrink the corrected area by a 1.05 factor to optimize the 
correction. 

To simulate the science image / (cf. Eq.Q} in polychromatic 
light, we sum 5 monochromatic images corresponding to wave- 
lengths uniformly distributed in the considered bandpass iS.A = 
Ao/Rj. Shannon sampling is imposed for fringes (/lo/£o)pix = 4, 
and £o follows Eq. [2] When photon noise is simulated, we con- 
sider a G2 star at lOpc observed by a space telescope with 40% 
of throughput and a bandwidth of 20% at Aq = 800 nm. We indi- 
cate the exposure time and the telescope diameter case by case. 
Table |2]i"egroups the other assumptions, which change from one 
section to another. Only one parameter generally takes realistic 
values so that its impact can be determined independently to the 
other parameters, set to ideal values. At the end of Sect. [9] all pa- 
rameters are set to realistic values to determine the SCC-FQPM 
performance. 



3.2. Criteria 

To optimize the different parameters of the self-coherent camera, 
we define two criteria. 



3.2.1. Averaged contrast 

The first criterion, called C\, gives the averaged contrast 
achieved in the corrected area <H of the SCC science image. This 
area varies as a function of three parameters : deformable mirror 
size, polychromatism and amplitude aberrations. Considering 
an N ac ,xN ac , deformable mirror: 

- in monochromatic light and without amplitude aberrations, 
<H is the centered (0.9/1.05 A/ acr )x(0.9/1.05 N act ) (A /D) 2 re- 
gion because we reduce the corrected area by shrinking it by 
a factor 1 .05 (Sect. 13. U and do not account for borders (fac- 
tor 0.9) of that area where the detection is prevented by the 
diffracted light of the uncorrected speckles. 

- in polychromatic light and without amplitude aberrations, 
following results of Sect. 18. II we reduce H to its intersection 
with the line of width and of same direction as fringes. 

If amplitude aberrations exist, we consider the region where the 
correction is effective (Sect. [7jl, which corresponds to half the ti 
area defined herebefore. The C\ can be written as 



Ci = 



I(a)6a 
/oJ K da ' 



(9) 



where Iq is the maximum intensity of the star image without 
coronagraph or SCC, and / is given by Eq.Q] The smaller C\ , the 
better the correction and the fainter companions can be detected. 
Since C\ is an average over the corrected area, it does not give 
any information on preferential positions in the image for the 
detection. 

As explained in Sect. 12.2.21 several iterations are required 
to reach very high-contrast imaging. That is why, when C\ is 
considered, we study its convergence speed as a function of the 
number of correction steps. 

3.2.2. 5 o- detection 

The second quantity we use to compare the different configura- 
tions is the 5 cr detection d$ a , which we define as 



5cr(p) 



(10) 



where cr(p) is the azimuthal standard deviation of the consid- 
ered image at radial separation p. As for C\, the smaller d^a-, 
the better the correction and the fainter companions can be de- 
tected. The dscr metric is a function of the angular separation in 
the considered image. It is then possible to determine preferen- 
tial angular separations for the detection. However, d^o-ip) does 
not give any information on favored directions. If dsAp) is used, 
we plot it against the angular separation. 

4. Intrinsic deformable mirror limitation 

To derive SCC performance, we needed to know the intrinsic 
limitation of deformable mirrors that we use. In that section, 
we set out that intrinsic limitation (no SCC) if no apodization 
is used. We considered a single telescope associated with a per- 
fect coronagraph and a deformable mirror. All the assumptions 
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section 



parameter 



m 



19.3.1.11 



19.3.2.11 



19.3.2.21 



SCC 



yes 



yes 



yes 



yes 



yes 



yes 



yes 



yes 



20 



1000 



1000 



1000 



1000 



10 and 20 



20 



25 



40 



16 to i 



16 to i 



150 



16, 32, 64 



32 



32 



64 



64 



16, 32, 64 



32 



32 



64 



coronagraph 



perfect 



perfect 



perfect 



perfect 



perfect 



perfect 



FQPM 



FQPM 



FQPM 



<r((/>) rms (nm) 



20 



20 



20 



20 



20 



20 



20 and 40 



20 



20 



tr(a) rms (%) 



Oto 1 
(T 3 PSD) 



1 

(fiat PSD) 



Photon noise 



yes 



Table 2. Simulation assumptions for each section. Notations are defined in Table 



7] Details as photon number are given in the text. 



are recalled in Sect. 13. II We simulated the coronagraphic image, 
without SCC, with a full DM correction (perfect estimation of 
the phase errors). We ploted d^ against the angular separation p 
in Fig.[3]for three different DM sizes and the same phase screen. 
The residue after the perfect coronagraph (no DM) is overplot- 
ted. 




1 i i i i I i i i i I i i i i I i i i i I i i i i I i i i i I i i i i 

5 10 15 20 25 30 35 



Angular separation in X/D 

Fig. 3. 5 o~ detection versus angular separation in the image pro- 
vided by a perfect coronagraph associated to a single telescope 
in monochromatic light without DM (solid black line) or with 
DM: 16x16, 32x32 and 64x64 from top to bottom. The initial 
phase error rms amplitude is 20 nm with a power spectral den- 
sity in /~ 3 . Residual phase error levels are given for each DM 
in nm. 



As expected, the cutoff frequency p c equals (N ac ,/2) (Aq/D), 
and higher frequencies are not corrected: no change for dsa-(p) 
with or without correction further than p c . For p < p c , the detec- 
tion limit depends on the diffracted energy of uncorrected speck- 
les located outside the corrected area (p > p c ). The larger N act , 
the larger pc and the lower the level of uncorrected speckles. 



That is why, for a fixed p < pc, the correction is better with a 
64x64 DM than with a 32x32. 

To improve the detection limit we can increase the number 
of actuators (degrees of freedom of the correction). However, 
nowadays, manufacturing issues limit that number to 64x64 
for MicroElectroMechanical Systems (MEMS). For the E-ELT 
planet finder, EPICS, a 200x200 device is under study but has 
not been developed yet. Another way to improve the detection 
limit is to minimize wavefront errors but it would be very opti- 
mistic to assume less than 20 nm rms for telescopes larger than 
a few meters. Another solution to rule out the detection limit is 
apodize the pupil so that the unco rrected speckles wou ld spread 
less light into the corrected area. iGive'on etalJ (2006) assume 
a Kaiser function to apodize the pupil and achieve a 10~ 1() con- 
trast (la detection) at 4Ao/D in the focal plane. But such an 
apodization pr ofile is are not easy to produce. W i th an am- 
plitude mask (IVanderbei et all 120031: iKasdin et all 120051) . its 
throughput would be about 15%, which significantly increases 
the exposure time for detecting faint companions. Another 
apod ization technique, c alled phase-induced amplitude apodiza- 
tion (iGuvon et alll20 05), could be used because its throughput is 
about 100%. But manufacturing p roblems limit i t s per formance, 
even if solutions are under study (Pluz hnik et al 1120061) . 

Finally, if no apodization is used, the best DM correction 
shown in Fig. [3] 

- cannot be overcome by any device used to estimate wave- 
front errors because it is an intrinsic limitation depending on 
the diffraction of uncorrected speckles; 

- comes from a speckle residue so that the sole way to improve 
it is a post-correction of the image by differential imaging 
techniques. One example is the self-coherent camera (step B, 
Sect. |23]l. 

We provide below results of a parametric study of the self- 
coherent camera used as a focal plane wavefront sensor. 

5. Signal-to-noise ratio 

We study the impact of the reference flux on the signal-to-noise 
ratio of wavefront error estimation (step A, Sect. 15. U and com- 
panion detection (step B, Sect. 15.2b . 

5. 1 . Step A: wavefront error measurement 

The purpose of step A is to measure wavefront errors directly 
from the noisy science image /. The interesting quantity is 
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then As and, here, we work out the signal-to-noise ratio of that 
measure considering the photon noise. 

To estimate As, we use Eq. [3] and obtain in monochromatic 
light 



A S ,n = 



(ID 



where the n index refers to the noisy quantities and where 7 is 
the Fourier transform of one of the two lateral correlation peaks 
of the science image /. Assuming A R is almost uniform in the 
image (see Sect.[6]l, EqfTTIimplies 



Var[A s ,„] 



(12) 



where we have replaced |Ar| 2 by the intensity of the reference 
channel on the detector 7 R . Considering photon noise, variance 
of 7_ n can be approximated by (see Eq.[42]of Appendix A): 



Var[I. in ] - 



(13) 



where (Ao/D) p i x is the number of pixels per resolution ele- 
ment Aq/D in the science image. From Eqs.[T2]and[T3] if the As, n 
measure is done with a signal-to-noise ratio higher than 77, we 
can write 



I/-., 



/r /^o\ 
1+/r \d) p 



>TJ. 



(14) 



Supplanting |7_| by |A S A R | = VTsTr (Eq.QI} and |A R | by a/Tr, 
we find: 



V/ s +/ R lD/ P ix ' 



(15) 



We want to estimate for wavefront errors (directly linked to 7s) 
and we know the reference intensity 7r. In the case of an asso- 
ciation with a FQPM (Sect. 19. 3.1. lb , 7 S » 7 R before beginning 
the correction loop. Under that condition, 77 max ~ (Ao/D)p lx V/r- 
Using the Shannon fringe sampling ((/lo/£o)pix — 3) and the min- 
imum separation between pupils of image and reference chan- 
nels (Eq.[2]i, the signal-to-noise ratio at the beginning of the loop 
roughly equals 5 VTr- After a few correction steps, we hope to 
achieve 7 S ~ 7 R . In that case, ?7 max ~ (A Q /D) pix V7 R /2 ~ 3.5 a/Tr. 
Finally, wavefront estimation (step A) can be done with a good 
accuracy as soon as 7r is greater than a few photons per pixel. 
The sole problem would be a zero value of the reference flux 
because the speckles would not be encoded and 77 would be 
zero (see Sect. [6). 

5.2. Step B: companion detection 

Step B corresponds to the companion estimation 7c est from the 
science image 7 of the last iteration of the correction (end of 
step A). It would be quite exhausting to study the exact propa- 
gation of noises through all steps of the algorithm that we use to 
compute the estimator of Eq. [7] We may present such a study in 
a future paper but it is not the purpose of the present one. Here, 
we express the dependence of the signal-to-noise ratio of the es- 
timator on fluxes in the image channel (7s, speckle noise) and in 
the reference channel (7r). 

We first assume uncorrected noises for the unmodu- 
lated (7 centn ) and modulated (|7_ n | 2 ) parts of the science im- 
age 7 (not rigorously exact). Considering the recorded reference 



intensity 7r is almost uniform on the detector (see Sect. |6), we 
can write the variance of the estimator 7c es t,n (Eq.[7]i: 



Var[I c est.n ,nj ' 



Var[\IJ 2 ] 



(16) 



We restrict the covariance of 7_ n to its variance (see Eq. [40] of 
Appendix A) and we suppose the distribution of the noise on 7_ n 
still follows a Poissonian law. We can then write the fourth mo- 
mentum of 7 n : 



Var[|7_, n | 2 ] - Var[|/_. n |] + 3 Var[|/ >n |] 2 . 



(17) 



Finally, using Eqs. [42] and [43] of Appendix A (variances of 7_ n 
and 7 cent n ) and Eqs. [T6land[T7l we obtain an approximation of 
the noise on 7c es t,n: 



Var[I Cesl , n ] - 7(a) j£j 



pix 



1 37(or) IDY 



ll ^ t W 



pixy 



(18) 



where 7 roughly equals to 7s + 7r + 7c. We insist that this 
expression is not exact, but we did numerically check that it 
roughly approximates the real noise. We can then derive the 
ideal flux regime (7s and 7r) to optimize the speckle substraction 
when applying estimator of Eq. [7] To do it, we assume the 
companion intensity is low compared to the coronagraphic 
residue (7c <s 7s). If this is not the case, the companion is 
detected and the SCC not useful. 

Weak fluxes 

If 7c <K 7s — 7r <k 1 photon per pixel, the variance of the noise 
on 7cest,n reduces to 



2 I D\ 

^r[7 Ce st,n] ^ - - ' 

7r \^o/ pix 



(19) 



and the signal-to-noise ratio /3 of the companion detection is 



/3-7 c 




«c 1, 



(20) 



which is not interesting. 

Dominating coronagraphic residue 

If 7 C <K 7 S , 1 <<K 7 S and 7 R <sc 7 S , Eq.[T8lbecomes 

37 2 /m 4 

Var[7 C est,n] - — r \ -^J Z 4 



'R 



(21) 



' pix 



If 7r < 1 , noise on the companion estimator is greater than the 
initial speckle noise and the SCC is useless. If 7r > 1, speckle 
noise is lightly reduced and the SCC acts positively. 

Dominating reference image 

If 7c <s 7s <K 7r and 1 <K 7r, we simplify Eq.[T8las 



Var[7 Ces t,n] ^47 R a 7 R . 



(22) 



The reference-image photon noise dominates the estimator 
noise. As it is greater than the speckle residue that we want to 
reduce (7s <k 7r), this case is not interesting again. 
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Strong fluxes for both channels 

If Ic <s h and 7r ^ 1$ » 1 photon per pixel, the variance of the 
noise on the companion estimator 7c es t,n (Eq.[T8ll reduces to 



V<»1/i-,st.„l - 8 7s oc/s, 

^O/pix 



(23) 



which means that it is proportional to the photon noise of the 
coronagraphic residue (of which the variance is 7s + 7r 2 7s): 
the speckle noise is suppressed. The signal-to-noise ratio ft of 
the companion detection is then 



7c lh\ ^ 37c 
8 7 s \D / P ix ~ 7 S 



(24) 



We have assumed (/lo/7))pj x ^ 5 to verify the Shannon sampling 
of fringes ((/lo/£o)pix - 3) and the minimum separation between 
image and reference pupils (Eq. [2j. The factor 3 reveals that the 
noise is smoothed and averaged when we select the correlation 
peaks and work out 7_ „ and 7 cen t,n by fast Fourier transform. 



we would like to set the smallest 7)r that we could. However, 
in a real setup, the light entering the reference channel is ex- 
tracted from the coronagraph rejected light, and its flux has a 
finite value. A minimum value for 7)r is then required to verify 
the condition of Sect. 15.11 and a trade-off has to be derived for 
each setu p. In the case of a n association with a FQPM coro- 
nagraph (Rou an et all 12000). we find that a reasonable value 
of y = D/Dr is between ~ 10 to ~ 30 in function of wavefront 
error level (Sect. |973l , and we often use y — 20 in the paper. 

We now quantify the impact of errors on the recorded ref- 
erence amplitude Ar. We first consider tip-tilt variations in the 
reference channel that induce spatial drifts of the reference im- 
age on the detector between Ar and Ar. With y — 20 and un- 
der the assumptions mentioned in Sect. 13.11 we find that for 
a ~ 5Aq/D (respectively 6) drift of the reference image, the 
correction is effective and roughly converges to the DM limi- 
tation (twice the DM limitation). The specification is then not 
critical, and in Sect. [9] we propose a compact, robust, and very 
simple device where the reference is stable enough in time. 



Finally, the best way to estimate the companion image 
from Eq. [7] with a good signal-to-noise ratio is to set strong 
fluxes in both image and reference channels (7s and 7r larger 
than 1 photon per pixel on the detector). 

For both wavefront estimation (step A) and companion es- 
timation (step B), the range of accepted reference fluxes 7r is 
broad and we propose to work with 7r around a few tens of pho- 
tons per pixel in the science image. 

6. Reference beam 

Like all differential imaging techniques, SCC needs the record- 
ing of a reference whose the complex amplitude Ar is used to 
estimate wavefront errors (Eq. [6j>. A first problem (zero division) 
appears at the positions where Ar is zero because the condition 
of Sect. 15.11 is not verified and the corresponding speckles are 
not encoded well. The second problem is that it is impossible 
to simultaneously measure Ar and the science image. We then 
propose to record the reference complex amplitude before be- 
ginning the loop and we call it Ar. But in that case we need that 
recording to be stable in time. In Sect. 16.11 we present a way 
to answer both problems at the same time. Then, we study the 
impact of errors on Ar: spatial drift 16. II and optical path differ- 
ence variation (Sect. 16.2b between the recording of Ar and the 
beginning of the loop. 

It is important to distinguish the reference amplitude Ar used 
to simulate the SCC science fringed image (see Sect. 13.11 ) and 
the reference amplitude Ar used in the estimator of Eq. [6] The 
second one is an estimation of the a priori unknown first one. 
In this paper, we use the diffracted complex amplitude by a 7)r- 
diameter pupil free from any aberrations for Ar. Flux and posi- 
tion of Ar are set using a recorded image of the reference chan- 
nel before the beginning of the correction loop. 

6.1. Reference pupil diameter and spatial drifts of the image 

The solution proposed to stabilize the reference channel and 
avoid the zero divisions is to use a small diameter 7)r for the 
reference pupil. In that way, Ar, roughly equal to the central 
part of the Airy complex amplitude, has very low sensitivity to 
wavefront variations in the reference channel, and its first dark 
ring is pushed away from the center of the image. At this point, 



6.2. Optical path difference 

Another potential error on the estimated reference amplitude Ar 
is the variation of the optical path difference (OPD) between the 
image and the reference channels. If the OPD varies, the cen- 
ter of the fringe pattern is shifted on the detector and, if we do 
not account for it, the wavefront error estimation is degraded. 
To quantify that limitation, we consider that the recorded ref- 
erence Ar corresponds to a zero OPD and we simulate cor- 
rection loops for several non-zero OPD under assumptions of 
Sect. I3.ll Figure |4] gives the convergence of the averaged con- 
trast C\ in the corrected area. If the OPD is zero, the correction 



monochromatic, 7=1000, DM32x32 




10~'° l i ■ ... i ■■■■■■■■ ■ 

5 10 15 

number of correction steps 

Fig. 4. Averaged contrast C\ in the corrected area against cor- 
rection steps for different optical path differences (OPD) be- 
tween reference and image channels. We assume a 32x32 de- 
formable mirror. 
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reaches the DM limitation (Sect. |4j) in 3 steps. When increas- 
ing the OPD, the correction converges more and more slowly. If 
the OPD is around the critical value /lo/4, the correction does not 
converge anymore. In that case, the estimated reference ampli- 
tude Ar (used in the wavefront estimator) equals to the real ref- 
erence amplitude Ar multiplied by i (imaginary unit). All phase 
errors are estimated as amplitude errors and vice versa. The cor- 
rection is not effective. To have the correction converging close 
to the DM limitation, we need to know the OPD between ref- 
erence and image channels with an accuracy of ~ /lo/6. And 
we have to stabilize it with the same precision (~ Aq/6) dur- 
ing the correction loop. It is possible to determine the OPD with 
accuracy just before closing the loop by recording a science im- 
age with a large bandwidth. The white fringe corresponds to the 
zero optical path difference, and its spatial shift from the opti- 
cal axis is directly linked to the OPD. To ensure the stability, we 
propose a robust compact device in the case of the association 
of the SCC with a coronagraph using a Lyot stop as the FQPM 
coronagraph (Sect. |9T3.1.21 >. 

7. Amplitude aberrations 

In that section, we look at the impact of amplitude errors in 
the pupil plane upstream of the coronagraph. We simulated the 
correction loop for several amplitude aberration levels (0% to 
1% rms). With a sole DM, amplitude aberrations induce a re- 
duction of the corrected area by a factor 2 and speckles of the 
uncorrected half-area diffract their light into the corrected half- 
area as seen in Fig. [5]- 1% amplitude aberrations and other as- 
sumptions detailed in Table [2] In Fig. [6] we plot the 5 cr detec- 




64X/D 

Fig. 5. SCC Science image after 4 steps of correction. Phase er- 
rors are set to 20 nm rms and amplitude aberrations to 1% rms 
before correction. 

tion d$ u- (Sect. 13. 2t at the 4 th iteration versus the angular sepa- 
ration for several amplitude aberration levels - see assumptions 
in Table [2] The best correction depends roughly quadratically to 
the amplitude aberration level. For example, with 1 % rms ampli- 



Amplitude aberration impact 




5 10 15 20 25 30 35 

Angular separation in A /D 



Fig. 6. Detection at 5 cr dso- against the angular separation in 
the image at the 4 th iteration of the correcion loop. We recall the 
coronagraphic residue without any correction (top black solid 
line). 

tude errors, the 5 cr detection is limited to ~ 1CT 7 at ~ 5 Ao/D. It 
is quite important to notice that this effect is a general limitation 
independent from the self-coherent camera. Whatever the tech- 
nique used to estimate for the phase and amplitude errors, it will 
be required to control with high accuracy the amplitude aberra- 
tions - less than ~ 1/1000 rms seems to be a minimum with 
an /~ 3 PSD - to achieve the performance limited by the DM. 
As explained in Sect. [4] using an apodized pupil might relaxes 
that specification. However, apodizing a pupil introduces other 
issues (throughput or manufacturing limitation) and a trade-off 
is needed. 

8. Chromatism impact 

The point spread function (PSF) size and the interfringe are both 
proportional to wavelength (cf. Eq.[TJ. In white light (small = 
/lo/AT), the science image is the superposition of all the 
monochromatic images over the considered bandwidth and 
fringes are blurred as seen in Fig. [7] As a consequence, speck- 
les are not encoded well and cannot be corrected with accuracy. 
Here, we present how chromatism limits wavefront estimation 
and we propose a hardware solution to use the self-coherent 
camera in polychromatic light (Rj = Aq/AA ^ 5 - 10). In 
Sect. 18.11 we study the field limitation for large bandwidths. In 
Sect. 18.21 we empirically optimize the chromatic factor F ap- 
plied in Eq. [6] Then, in Sect. 18.31 we quantify how the aver- 
aged contrast in the corrected area is modified when the band- 
width is enhanced. Finally, we describe a hardware solution to 
use the SCC in large bandpass in Sect. 18.41 

8.1. Field limitation 

We call A max and A m j n the largest and the smallest wave- 
lengths of the spectral bandwidth Ra- The mean interfringe 
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Fig. 7. Science images in monochromatic (left) and polychro- 
matic (Ra — 10, right) lights. These images correspond to the 
iteration (no wavefront correction). 

is /lo/^o (see Eq. |2). We assume the white fringe - null OPD 
- is in the center of the image. We derive the distance where 
the fringe systems for A max and A m j„ are shifted by half an inter- 
fringe: 



with n = as £o/Amax- We deduce n = Ra/2 and for Ra > 1: 

2 £o 



(26) 



In the perpendicular fringe direction, fringes become blurred 
as from a%. In the fringe direction, fringes are not blurred 
since OPD between the two channels is the same for all wave- 
lengths of the bandwidth. In that way, speckles of the science 
image are correctly encoded in the fringe direction whereas they 
are not in the perpendicular fringe direction as from a%. This 
effect is visible in the Fig. [8] where we show the 10th itera- 
tion of the correction for a 32x32 deformable mirror and two 
bandwidths. For these images, we do not apply the chromatic 
factor F(a) (Eq. |28| i in the wavefront estimator (Eq. [6]). Even 
if fringes are not blurred in the fringe direction, the speckles 
are dispersed (PSF size is wavelength dependent), and wave- 
front errors are not estimated perfectly. Hereafter, we study 
the combined impacts of both effects (fringe and speckle dis- 
persions). Chromatism limits the corrected area as soon as 




Fig. 8. Science image for the 10th iteration of the correction for 
two bandwidths: Ra = 20 (left) and 100 (right). We assume 
a 32x32 deformable mirror. The fringe direction is from bottom 
left to top right. 



is smaller than the DM cutoff spatial frequency, in the worst 
case V2pc with V2 there to account for the fringe direction). 



Using an A^ act xA^ act deformable mirror, pc = (iVact/2) Ao/D. We 
deduce from Eq. [26] the minimum spectral resolution RAmm to 
have no blurred speckles in the corrected area: 



R 



A min 



JV ac 



V2^. 
D 



(27) 



where £ ( > is the separation between the reference and the im- 
age pupils (cf. Table [TJ. If we consider^) - 1.65D (Eq.|2]i, we 
find R^mm — 37, 75, and 150, respectively, for a 16x16, 32x32, 
and 64x64 deformable mirror. 



8.2. Chromatic factor 

The wavefront error estimator given by Eq. [6] requires knowl- 
edge of th e chromatic factor F that we optimize in the current 
section. In lGalicher et al.l d2008l) . we used F = F t h derived from 
the model of light propagation through SCC: 



Fihi 



a) - f -^jexp(2ina.^Q/A)dA. 
Jr. A 



(28) 



However, F± (close to cardinal sin function) has zero values in 
the corrected area as soon as ag < pc (^th(«B) = 0). Speckles 
located at the zeros of F& are not encoded well and can be nei- 
ther estimated nor corrected. This limitation is similar to the zero 
values of the reference amplitude (Sect. 16. il l. Using the estimator 
of Eq. [6] pupil spatial frequencies corresponding to separations 
higher than (blurred speckles) cannot be recovered. But we 
can restrain their impact on the corrected area, and we test five 
chromatic factors: 



(29) 





= F a, 
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(F th + Fl) > 



We attempt to correct for unblurred speckles (Sect. 18. il l and min- 
imize the impact of the uncorrected speckles (blurred ones). To 
compare the different chromatic factors, we consider assump- 
tions given in Table [2] and look at the evolution of the aver- 
aged contrast in the corrected area C\ during the correction 
for Ra - 40. We find that setting the factor F to 1 /F t h provides 
the best results (no divergence). That factor reduces the impact of 
the uncorrected speckles because it multiplies their intensity by 
a low value in the estimator of Eq. [6] It acts as a regularization 
of the uncorrected spatial frequencies. We use that chromatic 
factor in all the paper from now. A more sophisticated estima- 
tor (regularization of a x 2 minimization without linearization of 
wavefront errors and without assumptions on Ra) could certainly 
be developed but has not been studied yet. 

8.3. Correction level 

The convergence speed of the correction loop slightly decreases 
when the light becomes more and more chromatic: 3 steps in 
monochromatic light and ~ 10 steps for Ra — 30 for a 32x32. 
At the same time, the averaged correction gets worst as seen 
in Fig. [9] where the criterion C\ at the \0th loop iteration is plot- 
ted against the spectral bandwidth (Ra = 20 to oo) under the 
assumptions of Table [2] In quasi-monochromatic light, the DM 
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Fig. 9. Averaged contrast C\ in the corrected area at the 10th 
iteration of the correction versus the spectral resolution R^ for 
three deformable mirrors. 



limitation is achieved (cf. Sect.|4]i for all the DM sizes. Naturally, 
the correction is more efficient when the number of actuators 
increases. The minimum Rj required to have the correction con- 
verging to less than twice the DM limitation in fewer than 5 steps 
increases with the size of the corrected area: ~ 80 for the 16x16, 
~ 100 for the 32x32, and ~ 200 for the 64x64. These values are 
close to the expected ones determined in Sect. 18.11 

From these results, it is clear that chromatism strongly lim- 
its SCC performance and that the self-coherent camera cannot 
be used with a classical bandwidth Rj ^5. Working with a nar- 
row bandwidth is not the solution because of photon noise. We 
foresee two possibilities. We can develop software solutions and 
modify the wavefront estimator of Eq. [6] to account for poly- 
chromatic dispersions of speckles and fringes (regularization of 
&X 2 minimization without linearization of aberrations and with- 
out assumptions on R^). Hardware solutions are also conceiv- 
able. We could associate the self-coherent camera with an inte- 
gral field spectrometer (IFS) at modest resolution (R,\ ~ 100). 
We could estimate wavefront errors for each wavelength chan- 
nel. We could also develop a new algorithm to process the data 
of all the channels at the same time to optimize the estimation. 
This solution will be studied in future work. Another hardware 
solution is a Wynne compensator that we describe in Sect. 18.41 



8.4. Wynne compensator 

The image widths (about A/D and A/Dr) and the fringe pe- 
riod (A/go) are proportional to wavelength, which is the reason 
for the SCC chromatism limitat ion (Sec ts. [8T| to |8.3t . In t he con- 
text of speckle interferometry (LabeyriH ll970h 7 IWynnd d 19791) 
proposed a device to correct for such a spectral dependence over 
a wide spectral range: the Wynne compensator (Fig. [Tol l. This 
device is composed of two triplets of lenses. In each triplet, 
the first and the third lenses are made in the same glass, called 
glass 1, while the second lens glass is different, glass 2. Each 



Triplet lenses 
Science 



Triplet lenses 



beam 





Reference 
beam 

Fig. 10. Scheme of a Wynne Compensator composed by two 
triplet lenses made with two kinds of glasses (called 1 and 2). 
Dispersion is linear with the wavelength and the outgoing beam 
is collimated. 

second lens is symetric (same curvature on each side). Indexes 
of glasses 1 and 2 are equal at the middle wavelength of the con- 
sidered bandpass, but their dispersions are different. The beam is 
then not modified for the middle wavelength (black solid curve), 
while it is magnified for other wavelengths (red dotted and blue 
dashed curves). The compensator's outgoing beam is collimated 
for all wavelengths. In Fig.[l0] we only show the average beam 
of image and reference channels for reasons of clarity. In that 
way, we only outline that the Wynne compensator magnifies 
the separation t;o(A) between pupils. A more detailed study of 
the Wynne compensator shows it also magnifies pupil diame- 
ters D(A) and Dr(A). By choosing the right glasses and opti- 
mizing the curvatures of the lenses, their thicknesses and the 
distance between the two triplets, we can apply a magnification 
proportional to the wavelength over a large bandwidth in visible 
light. 

Such a Wynne compensator can be associated to the SCC by 
adding it just before the recombining optic. Since magnification 
is proportional to the wavelength, D, Dr, and are proportional 
to A and we obtain non-blurred fringes all over the detector (cor- 
responding to the monochromatic case). A simulation is shown 
in Fig.[TT]for a spectral resolution of Rj = 6.5 at 650 nm. Since 
only magnification is important, we do not consider any corona- 
graph and the incoming wavefront is assumed to be aberration- 
free. The relative difference between the Airy pattern widths and 




Fig. 11. Science image without (left) or with (right) Wynne com- 
pensator. The initial spectral resolution is 6.5 at 650 nm. For the 
corrected image, the effective bandwidth is about 0.75% ( right). 
We use log scale. 

interfringes for the extreme wavelengths (600 and 700 nm) is 
about 0.75%, which corresponds to Ra — 130. That spectral res- 
olution (130) is the spectral bandwidth required to have almost 
no blurred speckles in the corrected area of a 64x64 deformable 
mirror (Sects. [8TI and [Oi l. An optimized Wynne compensator 
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would enable us to use the self-coherent camera in polychro- 
matic light (R \ = 6.5 in visible light Aq = 650 nm) and achieve 
the quasi-monochromatic performance (Fig. [9]). 

Finally, even if chromatism seems to be a hard point of the 
self-coherent camera technique, several solutions are conceiv- 
able: more sophisticated estimators, hardware solutions (associ- 
ation with a Wynne compensator or an integral field spectrome- 
ter at modest resolution). 



9. SCC and real coronagraphs 

The previous sections present a parametric study of the self- 
coherent camera without mentioning a ny concrete se t up an d 
as suming a perfect co ron agraph as in iGalicher et al] ([2008 ). 
In lBaudoz etail d2006l) and lGalicher & Baudozl d2007l) . we pro- 
posed a self-coherent camera device built as an interferome- 
ter: a beamsplitter to create the reference channel, a pinhole 
to filter the reference beam, a delay line to ensure a null op- 
tical path difference, and a lens to recombine image and ref- 
erence beams. The disadvantage of that device is the delay 
line that has to be controlled with very high accuracy in real 
time (Sect. |6~2l ). In Sects. [9~T| and |9. 21 we propose a new robust 
design for associating the self-coherent camera with a corona- 
graph that h as a Lyot stop plan e. We describe in detail the case 
of a FQPM dRouan et aUl2000l) in Sect, f 



9.1. How to use the Lyot stop plane 

Using a coronagraph that needs a Lyot stop is very interesting 
because it can be easily associated with the SCC. Such a coron- 
agraph rejects only the stellar light outside the pupil - no com- 
panion light - so that the light stopped by the Lyot stop comes 
only from the hosting star and can be used to create the refer- 
ence channel. We propose adding a small non-centered pupil to 
the classical Lyot stop as shown in the schematics of Fig. Q~2] 
for a FQPM coronagraph. The sole difference from the classical 




Science 
pupil 



SCC Lyot stop 



Fig. 12. Schematics of the self-coherent camera associated to a 
coronagraph that uses a Lyot stop plane. The pupil intensity dis- 
tribution in the Lyot stop plane is given for a FQPM corona- 
graph. 

coronagraphic device is the modification of the Lyot stop. That 
solution is very attractive because it constrains a constant optical 



path difference between the reference and the image channels. 
As explained in Sect. 18.41 we may w ant to associat e the SCC 
with a chromatic Wynne compensator (lWvnneLll979l) . This one 
would be placed between the SCC Lyot stop and the recombin- 
ing optic. 

9.2. Wavefront estimation 

The complex amplitude estimated by Eq. [4] is the pupil ampli- 
tude in the Lyot stop plane downstream the coronagraph. The 
second step is determine the complex amplitude tf/' s upstream of 
the coronagraph where the deformable mirror stands. If the coro- 
nagraph is perfect and wavefront aberrations are small, we can 
use the estimator of Eq[6j If the coronagraph is not perfect, we 
have to model the light propagation through it. 

Calling M the mask function in the focal plane and L the 
classical Lyot stop (sole image channel), we have 



<As =(^*r" 1 (M)) L, 



(30) 



where * denotes the convolution. We deduce \f/' 5 within the L 
area: 



-As = r 



(31) 



where (l/M)o is the inverse of the mask function M where M 
is not zero and equals elsewhere. This expresses that we can- 
not estimate the spatial frequencies for which the mask M has 
stopped the energy (ie. M = 0). We finally estimate the complex 
amplitude upstream of the coronagraph from Eqs.[4]and[3T] We 
notice that phase masks are not bounded by this limitation since 
they do not block light. 

9.3. SCC and four-quadrant phase mask 

The FQPM coronagraph uses a Lyot stop and can be associ- 
ated with the self-coherent camera in a device, called SCC- 
FQPM. The FQP M has been d e scribed in d etail i n pre- 
vious papers dRouan eta D, 120001: iRiaud et all l200ll 120031: 
iBoccaletti et alll2004l) . We recall the coronagraphic mask M in- 
duces a n phase shift on two quadrants of a diagonal and no 
phase shift in the two other quadrants. In Fig. [12] we explain 
where the image and reference pupils are picked in the Lyot stop 
of an SCC-FQPM. In Sect. 19.3. 1.11 we first give the quantity of 
energy that goes through the reference pupil to determine the 
optimized ratio y = D/Dr of pupil diameters (cf. Sect. 16. It . We 
then study the optical path difference between reference and im- 
age channels (cf. Sect. 16. 2t in Sect. 19.3.1.21 Finally, we examine 
the performance of the SCC-FQPM under realistic conditions in 
Sect. 19.3.21 In the whole section, we assume a perfect achromatic 
FQPM: infinitely thin transitions and perfect phase shifts. 

9.3.1. Implementation 
9.3.1.1. Reference flux 

As explained in Sect. I6.ll the ratio between image and ref- 
erence pupil diameters y - D/Dr has to be optimized for 
each SCC design. In Table [3] we give the averaged intensi- 
ties of the image and reference PSFs in the corrected area of 
a 32x32 DM for different y and phase-error levels. All fluxes 
are expressed in photons per pixel. To establish the table, we 
consider the device of figure Q~2] working under assumptions of 
Table [2] We put the reference pupil behind one of the quad- 
rants of the FQPM (see Sect. I9.3.1.2K and we set the stellar flux 
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using assumptions of Sect. 13.11 a 20 s exposure time, and 4 m 
diameter telescope. Coronagraphic residue < Is > in the im- 



r 


10 


20 


o-[</>] (nm) 


20 


40 


20 


40 


< Ir > 


212 


193 


58 


53 


< h > 


45835 


176218 


45835 


176218 


</s> 


0.46 


0.11 


0.13 


0.03 



Table 3. Averaged fluxes of image (< Is >) and reference ( < 
I R >) channels in the corrected area of a 32x32 DM for differ- 
ent y (first raw) and phase error levels (second raw). We con- 
sider an achromatic FQPM. Fluxes are given in photons per 
pixel. See Table\2\for more informations on assumptions. 



age channel depends roughly quadratically on the phase error 
level as expected because the coronagraphic pupil is approxi- 
mately fts — {a + i(f>)P. In contrast, the rejected star energy 
remains slightly constant in the reference pupil (< /r >). We 
notice also that < Ir > roughly evolves as the inverse square 
of y. We propose setting y to ~ 20 so that the condition of a few 
photons per detector pixel coming in from the reference chan- 
nel (see Sect. 15. It can be respected with exposure times that are 
not too long (speckle lifetime). It is important to notice that this 
choice is not strict. To strictly set y, we have to know the lifetime 
of the speckles to be corrected. This would impose a maximum 
exposure time to record the interferential images / (aberrations 
have to be static to be estimated from /). From that maximum 
time, we could choose the apropriate y value to let the reference 
flux verify the condition of Sect. 15.11 To make that choice, we 
may notice again that higher values of y are interesting to ensure 
stability of the reference image and to limit the presence of low 
values of Ar in the corrected area (cf. Sect. 16. lb . 



employ a more complex model accounting f or the exact FQPM 
impact on the first 999 Zernike polynomials dNolll [19761) . Using 
that model, we plot the averaged contrast C\ in the corrected 
area of a SCC-FQPM versus the correction step for several spec- 
tral bandwidths (Fig. [T3l , Assumptions are given in Sect. 13.11 
and Table [2] The same remarks as in a perfect coronagraph 
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Fig. 13. Averaged contrast C\ in the corrected area of a SCC- 
FQPM versus the number of correction iterations for several 
bandwidths Ra — Aq/AA. The dot-and-dash line and the dashed 
line represent respectively the best DM correction and the 
monochromatic case with a perfect coronagraph. 



9.3.1.2. Optical path difference 

We have shown in Sect. I6.2l that we need to accurately control 
the OPD between reference and image channels (accuracy and 
stability of ~ A /6). In the case of an SCC-FQPM, if the refer- 
ence pupil is picked behind one of the n phase shift quadrants of 
the focal mask, the optical path difference in the Lyot stop plane 
in monochromatic light is n. Fringes are shifted by half an inter- 
fringe, and we have to multiply by -1 the estimator of Eq.FJT] If 
the reference pupil is picked behind another quadrant, the OPD 
is null and we can use directly the estimator of Eq. [3J] The di- 
rections to avoid are the transition direction of the FQPM for 
which the OPD is not just or n and above all is not uniform 
over the reference pupil so that it is not easy to account for it 
in the estimator. In a real experiment, we propose to pick the 
reference pupil behind one of the four quadrants avoiding the 
transition directions and to calibrate the sign of the estimator. 
We are confident that the OPD is stable in time in the compact 
device proposed in Fig.Q~2]and laboratory demonstrations are in 
progress to confirm it. 

9.3.2. Performances 

9.3.2.1 . Impact of chromatism 

If we directly apply the wavefront estimator of Eq. [31] the cor- 
rection loop becomes slower than in the perfect coronagraph 
case (Sect. |8.3T l because the FQPM model is not perfect. We then 



case (Sect. 18 . 3b are valid: convergence slows down and gets 
worst when the bandwidth gets wider (Rx decreases). However, 
the SCC-FQPM best correction is around five times less effec- 
tive than with a perfect coronagraph (dashed line) which reaches 
the DM limitation (dot-and-dash line, cf. Fig.0. This is because 
the FQPM does not well estimate all optical defects. For exam- 
ple, the astigmatism in the FQPM transition direction is very 
poorly estimated (study not shown in this paper). We work on 
that limitation to overcome it, and we also plan to use an other 
coronagraph mask without such effects as the annular groove 
phase mask dMawet et all 120051) . However, even with that limi- 
tation, the correction is effective and SCC-FQPM provides very 
high-contrast imaging as shown in the next section. 

9.3.2.2. SCC-FQPM detections 

In the last section, we check for the detection efficiency of 
SCC-FQPM under realistic assumptions detailed in Sect. 13.11 
and Table [2] for spatial observations. The variance of the whole 
phase error is 20 nm rms, and we set the initial astigmatism 
defects in the FQPM transition direction to 1 nm rms (lev- 
els before correction). We account for amplitude aberrations 
of 1% rms. We simulated twelve 2 10~ 10 earths, located sym- 
metrically with respect to the image center at 3, 6, 9, 12, 15, 
and 18 Aq/D. For the 4 m telescope that we consider (throughput 
of 40%), these separations correspond to 0.12 (1.24), 0.25 (2.48), 
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Contrast (xlO 10 ) 


Angular separation in A /D 
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18,0 


18,0 
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7,7 


20,0 


19,9 


10 


9,1 


36,5 


36,3 



Table 4. Comparison between the contrasts and angular separations measured in the SCC image and the simulated values. 



0.37 (3.71), 0.50 (4.95), 0.62 (6.19), and 0.74 (7.43) arcsec (AU). 
One 810~ 10 super-earth is located at 20A a ID (8.25 AU) and 
a 10~ 9 Jupiter-like planet sits at 36.5 A /D (15 AU). All spec- 
tra are flat and the hosting star is a G2 star at 10 pc. We closed 
the loop to correct for wavefront errors (step A, Sect. 12.21 ) dur- 
ing 7 h25 min, which corresponds to 6 iterations. After step A, 
we apply the companion estimator and obtain the image pre- 
sented in Fig. Q3] We consider photon noise but no read-out 
noise. The white circle of uncorrected speckles corresponds to 



tions limitation). Intensities of the super-earth in the fourth quad- 
rant (blue circle) and of the Jupiter on the top right (red circle be- 
yond the reference ring) are well above the surrounded speckle 
field. Table [4] gives the measured flux and position for each de- 
tected planet. Fluxes are slightly underestimated (around 10% 
accuracy). Positions are accurate at ~ 0.25 Ao/D, which equals 
to lOmas (0, 10 AU for a lOparsec star). 

Finally, planets as faint as earths are detectable by SCC- 
FQPM in a few hours from space under realistic assumptions. 




Fig. 14. Detections of earths (white circles), super-earth (blue 
circle) and Jupiter (red circle) with a SCC-FQPM downstream 
a Am space telescope and a 64x64 deformable mirror. The image 
is the result of a wavefront correction loop (step A, 7 h25min) 
and a companion estimation (step B). 



the low values of the reference amplitude for y — 25. Stellar 
speckles at these positions are not correctly encoded (estimated) 
and they are not corrected (Sect. |5}. The red dashed box bounds 
the centered 64x64 (Ao/D) 2 area. Five earths are clearly de- 
tected (white thin circles) in the first quadrant at 0.25 (6), 
0.37 (9), 0.50 (12), 0.62 (15), and 0.74 (18) arcsec (A /D). 
The sixth at 3 Aq/D is detected (white thin circle in the fourth 
quadrant), but speckle noise is just below its level. Positions 
of the six other earths are symmetric with respect to the im- 
age center and are in the uncorrected area (amplitude aberra- 



10. Conclusion 

In Sect. H] we provided the intrinsic li mitation for deformabl e 
mirrors controlled by the algorithm of Borde & Traubl (|2006), 
under realistic yet optimistic assumptions (no dead actuators, 
continuous face sheet). It is important to keep in mind that 
this limitation does not depend on the technique used to esti- 
mate for wavefront errors since we assumed a perfect estima- 
tion. One way to improve the DM best contrast could be an 
apodization of the pupil so that the uncorrected speckles would 
diffract their light in a more restricted area. However, all the 
techniques p r opose d to apodize a pup il ( Vanderbei et al, 2003; 
iKasdin et al.L l2005t iGuvon et al.L l2005t Pluzhnik et ail [2006) 
come with throughput problems or manufacturing limitations. 

In Sects. |5]to[8] we gave the results of the parametric study 
of a self-coherent camera (SCC) used as a focal plane wave- 
front sensor and associated with a perfect coronagraph and a 
deformable mirror. Several points do not seem to be critical 
for the technique: reference flux (Sect. |5), error on the exact 
position of the reference image, and diameter of the reference 
pupil (Sect. 16. U . On the contrary, two points are more critical: 

- optical path difference between the reference and the im- 
age channels. We have to know and control this OPD with 
an accuracy of about Ao/6. If we associate the SCC with a 
coronagraph using a Lyot stop as described in Sect. [9] the 
hardwar optical path difference is always zero and we only 
have to control the end of the setup. Only common optics 
are used in this setup. We could put the device in a closed 
box to avoid differential air variations. We plan to check for 
the level of the variations in the optical path difference in the 
device of Fig.Q~2]in a laboratory experiment. 

- chromatism. As shown in Sect. [8] chromatism is the most 
critical point in the technique. The main consequence is the 
reduction of the corrected area, in other words, the field of 
view of the image as shown in Fig. [8] The uncorrected speck- 
les spread light and limit the contrast of the detection. We 
are studying software and hardware solutions to minimize 
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that effect. For example, for the former we could develop 
a more sophisticated wavefront estimator using a x 2 mini- 
mization with regularization terms to account for the spectral 
dispersion of the speckles (Aq/D) and of the fringes (£q/D). 
Hardware solutions are certainly more appropriate. We pre- 
sented the Wynne compensator in Sect. 18.41 According to 
numerical simulations, it would enable working with a clas- 
sical bandpass (~ 15%, Ra — 6) in the visible light with 
a 64x64 DM. We will test such a Wynne compensator in a 
laboratory experiment. A second hardware solution to over- 
come the chromatism limitation could be the association of 
the SCC with an integral field spectrometer at modest spec- 
tral resolution (R^ =2 100 - 150). This solution is very attrac- 
tive because it would directly provide a companion spectra, 
but no work has been done on it yet. 

We showed in Sect. [7] that the self-coherent camera can es- 
timate for both phase and amplitude aberrations. The impact of 
the amplitude aberrations is quite critical: their level must be 
smaller than 1 /1000 to reach a 5 <x detection of 10~ 9 at 5 Aq/D 
with an f~ 3 power spectral density. That limitation is not in- 
trinsic to the SCC and would limit any high-contrast imaging 
system. 

Section|9]presented a very simple and robust design that as- 
sociates the self-coherent camera with any coronagraph which 
uses a Lyot stop. The most interesting point is that the optical 
path difference between the two channels is constant per con- 
struction. In Sect. 19.31 we studied in detail the case of the as- 
sociation of the SCC with a FQPM coronagraph (SCC-FQPM) 
and we showed in Sect. 19.3.21 that the performance is very attrac- 
tive and comp arable to the case o f a perfect coronagraph that is 
presented in Galic heTet alj (12008). Detections of 2 10~ ll) earths, 
8 10~ 10 super-earth, and 10~ 9 Jupiter under realistic assumptions 
are numerically demonstrated for an SCC-FQPM in polychro- 
matic light (Ra = 5) using a Wynne compensator (reducing the 
effective bandwidth to /^eff = 150) in ~ 7h25min with a 4m 
space telescope. 

The next steps are laboratory demonstrations of both SCC 
capabilities: focal-plane wavefront estimation (step A, Sect. |2.21 > 
and companion estimation by differential imaging (step B, 
Sect. 12.3b . We will also attempt to overcome the poor estimate 
of the astigmatism in the FQPM transition direction. New algo- 
rithms have already been developed for using a DM interaction 
matrix including the impact of the whole instrument: DM, coro- 
nagraph, and SCC. It will be tested in a laboratory experiment 
very soon. We also study the SCC associati on with other coron - 
agraphs like an annular groove phase mask (Ma wet et aU l2005). 



We thank Remi Soummer for private communications 
about his pap er "Fast computation of Lyot-style coronagraph 
propagation" dSoummer et aUl2007bl) . which was very useful for 
simulating the polychromatic images and the different values of 
the y parameter. 

Appendix A 

In that appendix, we present how photon noise and read-out 
noise propagate through the numerical algorithms providing 
the wavefront estimation (Eq. |6j and the companion estima- 
tion (Sect. |521). We call I„(a) the noisy intensity of the recorded 
sience image: 



where 1(a) is the noiseless intensity at position a of the in- 
terferential image (Eq. [TJ and e(a) the noise. For photon 
noise (Poissonian distribution), the mean of e is zero and its vari- 
ance Var[e(a)] at the position a is 1(a). For a uniform read-out 
noise (gaussian distribution), the mean of e is null and its vari- 
ance is cr^ 0N over the whole image. For both photon noise and 
read-out noise, we assume the noise is not correlated between 
pixels located at different positions a\ and a 2 (the spatial covari- 
ance of e is zero): 



E[e(ai)e*(a 2 )] = Var[e(ai)]6(a l - a 2 ), 



(33) 



where 6 is the Kronecker symbol and E the statistical mean. 
The self-coherent camera technique uses < F~ l [I„], the inverse 
Fourier transform of the science image /„. Let v(u) be the noise 
on T~ x [I„](u) at the position u in the Fourier plane, such that: 



T-Vn\(u) = T-V](u) + v(u). 



(34) 



Using Eqs. [32] and [34] the monochromatic case at Aq and the 
Fourier transform properties, we find that the mean of v is zero 
and its spatial covariance is 



E[v(m) v*(u 2 )] = 

2in(a\.u\ 



ff 



exp 



a 2 .u 2 ) 



E[e(ai)e(a2)]da\da2- (35) 



From Eq. [33] we determine the covariance of v in the case of 
photon noise 



E[v( Ul )v\u 2 )] =T~ X [7](«i -u 2 ), 
and in the case of read out noise: 
E[v(ui)v*(u 2 )} = cr^ 0N 5(«i - u 2 ). 



(36) 



(37) 



The spatial covariance of the inverse Fourier transform of pho- 
ton noise (Poissonian distribution) is not reduced to its variance 
whereas it is for read-out noise (Gaussian distribution). 

During the SCC image processing, to estimate As, we iso- 
late one of the two lateral peaks of T x [I n ] dBaudoz et all 120061 
iGalicher & Baudozl 120071) using a circular binary mask of di- 
ameter D + Dr =2 D (see Sect. 16. II for that approximation) and 
we apply a Fourier transform. The noiseless part of that Fourier 
transform is called I-(a) (Eq.0). The noisy part is I-„(a) 



/_„(<*) = Ua) + e_(ar), 



(38) 



with e-(a) the noise at the a position. The mean of e-(a) is zero 
because E[v(u)] = 0. From Eq.[36] we find that in the case of a 
photon noise: 



£[e_(oi) el(a 2 )l 

r r r- l v]( Ul 

Jd Jd Aq 



112) 



exp 



2 in(a 2 .u 2 — a\.u\) 



A* 



d«i d«2 



,(39) 



where I represents the integral over a pupil of diameter D. 
Jd 

Calling yfo(a') = J\(a)/a where J\ is the Bessel function of 
the first kind corresponding to a D-diameter circular aperture, 
we can write 



I„(a) = 1(a) + e(a), 



(32) 



E[e-(ai)e*_(a 2 )] = 

jj I(a 3 )Mj> (Hon - orall) ^0 (llorj - a 3 ||) da 3 . 



(40) 
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We deduce from Eq.[40]the variance of e_ in the case of photon 
noise: 

Var[eja)] = Var[I^„(a)] = 1(a) * |^(||ff||)| 2 . (41) 

Practically speaking, we record the interferential image / on a fi- 
nite number of pixels and numerically process data. This means 
that the Fourier transform of Eq. [39] is a fast Fourier transform 
and the width of |^[| 2 is ~ (Ao/D) p [ x pixels. Finally, the |^l| 2 con- 
volution in the expression of Var[e_] roughly corresponds to av- 
eraging the noise over ~ (Aq / £>) 2 ix pixels and 

Var[/_,„(a)] - (42) 

A similar result is found for the noise on the unmodulated 
part 7 cen t of the image / Is + Ir + Ic (Eq.[8]i: 

4 1(a) 

Var[I cenUn (a)] - - ' ■ (43) 
(*o/D) 2 pix 

The sole difference is the size of the selecting binary mask in the 
correlation plane. 
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